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Abstract 

We continue numerical studies of the spectrum of the pure U(l) lattice gauge theory 
in the confinement phase, initiated in [1], Using the extended Wilson action S = 
— \j3 cos(0p) + 7 cos(20p)] we address the question of universality of the phase 

transition line in the (/3, 7 ) plane between the confinement and the Coulomb phases. 
Our present results at 7 = —0.5 for the gauge-ball spectrum are fully consistent with 
the previous results obtained at 7 = —0.2. Again, two different correlation length 
exponents, v ng = 0.35(3) and v g = 0.49(7), are obtained in different channels. 
We also confirm the stability of the values of these exponents with respect to the 
variation of the distance from the critical point at which they are determined. These 
results further demonstrate universal critical behaviour of the model at least up to 
correlation lengths of 4 lattice spacings when the phase transition is approached in 
some interval at 7 < — 0 . 2 . 
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1 Introduction 


The properties of the strongly coupled pure U(l) lattice gauge theory (U(l) 
theory) remain insufficiently understood after 20 years of its numerical investi¬ 
gation. Because of the two-state signal present on finite toroidal lattices (cubic 
lattices with periodic boundary conditions), it remains unclear whether a con¬ 
tinuum field theory can be constructed at the confinement-Coulomb phase 
transition. A two-state signal is obviously consistent with the first order, but 
it does not exclude the second order as long as the latent heat decreases with 
increasing volume and can be extrapolated to zero in the thermodynamic limit 
using a power law. Results of recent investigations [2-4] allow such an extrapo¬ 
lation in some range of coupling parameters. Furthermore, when homogeneous 
spherical lattices (as opposed to the surface of a 5-dimensional hypercube) are 
used, the two-state signal is not observed at all [5,6]. 

In such a situation there is little hope of answering the question of the order 
beyond any doubt while no reliable analytic means are available and lattice 
sizes in numerical simulations are restricted. Furthermore, the lattice action 
can be considerably varied and the extended Wilson action used in most works, 
and also in this paper, is only an example. There may exist regions of the pa¬ 
rameter space where the transition is clearly continuous. Indeed, in a recent 
work [7] the existence of a continuous phase transition possibly with an essen¬ 
tial singularity is indicated by results using a novel formulation of the pure 
compact Abelian gauge theory in continuum. 

Therefore it is appropriate to ask how far the properties of the phase transition 
are also consistent with the scaling behaviour associated with the continuous 
transition. If for a certain action such a consistency is found in some finite 
range of scales (in lattice cutoff units), two interesting scenarios are available: 
First, this consistency indicates the existence of a critical point not too distant 
in parameter space from the action used so that it dominates the behaviour 
of the chosen system in the observed range of scales. Second, if the range of 
scales is very large, the theory might be considered as an effective theory, with 
finite but large cutoff, even if no critical point is found. 

This may be considered a realistic approach to the question of whether the 
U(l) theory has a physical meaning because considerable indications of a 
second-order scaling behaviour for the extended Wilson action (parameterized 
by two simplest plaquette couplings f3 and 7) already exist. On homogeneous 
spherical lattices, a volume dependence of several bulk observables described 
by simple second-order finite size scaling theory was found in a large range of 
lattice sizes [5,6] at three points of the phase transition, 7 = 0, —0.2, —0.5, of 
the extended Wilson action and also for the Villain action [8]. This fact (and 
not the absence of the two-state signal on spherical lattices) is up to now the 
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most significant evidence for the existence of a non-Gaussian fixed point and 
universality in the U(l) theory. 


The other indication of second-order scaling behaviour has been found in the 
spectrum studies of the U(l) theory in the confinement phase on toroidal 
lattices. Working with the extended Wilson action at 7 = —0.2, the spectrum 
of the gauge-balls (GB) at various distances r = \j3 — (3 C \ from the phase 
transition was determined [1]. Its r-dependence has been found to be described 
in a very broad interval (outside the narrow region where the two-state signal 
appears) by simple power laws. Such a behaviour is predicted by the leading 
term in a second-order scaling picture, parameterized by the correlation length 
exponent v. Two different values of the exponent v were found, one of them 
non-Gaussian. This suggests the existence of an interesting fixed point in the 
pure U(l) theory. Of course, another possible explanation discussed in [ 1 ] is 
the occurrence of some cross-over or pre-asymptotic phenomenon. A possible 
way to distinguish between these possibilities is to look for universal properties 
and to repeat the calculations at a different point on the phase boundary. 


Therefore in this paper we perform a similar analysis at 7 = —0.5. The main 
aim is to check whether the results for v are consistent with those at 7 = 
—0.2. This turns out to be the case, further supporting the presence of a 
non-Gaussian fixed point, which is the most natural explanation of such a 
universal scaling behaviour. The alternative explanation that the results of 
[ 1 ] are influenced by a crossover related to a possible tricritical point in the 
vicinity is no longer probable. 


We pay particular attention to the stability of the values of v with respect to 
the variation of the /3-interval in which they are determined by power-law fits 
to the measured masses. We demonstrate, also for data obtained previously [1] 
at 7 = — 0 . 2 , that, within the lattice sizes and correlation lengths we could 
achieve, there is no significant indication that v approaches the value 1/4, 
expected for a weak first order phase transition. 


The new results are presented as follows: In section 2 we briefly describe the 
model and summarize our earlier results on the GB spectrum at 7 = —0.2. 
In section 3 we present results for the scaling behaviour of the GB masses at 
7 = —0.5. In section 4 the stability of the results at both 7 values is discussed. 
In section 5 we compare the results at both 7 values and draw our conclusions. 
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2 The U(l) theory and earlier results for GB masses at 7 = —0.2 


2.1 Action and phase diagram 


We consider pure U(l) lattice gauge theory with an extended Wilson action 


S = — ^2 [/?cos(@p) + 7 cos( 2 @p)]. ( 1 ) 

p 


Here, ©p e [0, 2n) is the plaquette angle, i.e. the argument of the product 
of U(l) link variables around a plaquette P, and (3 and 7 are the single and 
the double charge representation couplings, respectively. Taking ©p = a 2 gF fJjU , 
where a is the lattice spacing and P + A'y = 1/g 2 , one obtains for weak coupling 
g the usual continuum action S = \ / d 4 xF 2 u . 

This lattice gauge theory has a line of phase transitions between the strong 
coupling confinement phase and the weak coupling Coulomb phase. Its position 
for the Wilson action (7 = 0) is /3 C — 1.011 [9], whereas for 7 = —0.5 it is 
given below in (12). For 7 > +0.2 the transition is most probably of 1 st order, 
weakening with decreasing 7 [10,9]. 

The recent studies on spherical lattices suggest that the order changes at 
7 = 70 ~ 0, and is of 2 nd order for 7 < 70 [5,6]. With decreasing 7 the 2 nd order 
transition further weakens in the sense that the specific heat peak decreases 
for fixed lattice size and the autocorrelation time increases [ 6 ]. The scaling 
behaviour of bulk quantities is universal at least in the range —0.5 < 7 < 0. 

On toroidal lattices the disturbing two-state signal weakens with decreasing 
7 , but is present at least until 7 = —0.5 [9]. For even smaller 7 the large 
autocorrelation time makes simulations prohibitively expensive. Thus on the 
toroidal lattices the two-state signal on the critical line cannot be avoided. 


2.2 Summary of the earlier results at 7 = — 0.2 


Now we briefly summarize those results obtained for 7 = —0.2 in [1] which 
are relevant for this work. In the vicinity of the critical point 

(3 C = 1.1609(2) (2) 


in the confinement phase two groups of GB masses have been found with 
distinctly different scaling behaviour when the phase transition is approached, 
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r —> 0. Most of the GB masses scale proportionally to T Uns , the value of the 
correlation length exponent z/ ng being non-Gaussian, 

u ng = 0.35(3). (3) 


This agrees with the exponent v = 0.365(8) found in FSS studies on spherical 
lattices [5,6,8]. 

The mass of the J pc = 0 ++ GB scales as r^ g , with 

v, = 0.49(7). (4) 


This correlation length exponent agrees with the Gaussian value v = 1/2. We 
also give the result for the ratio of both exponents, 

= 0.71(8). (5) 


We note that an observation of two different critical exponents at one value 
of 7 leaves open the possibility that one of these exponent values is caused 
by crossover phenomena associated with the possible presence of a tricritical 
point in the vicinity. Two such scenarios have been discussed in [1], Their con¬ 
firmation or rejection requires repeating the calculations either at the same 7 
and substantially larger lattices, which is currently not possible, or at another, 
distant point in the parameter space. 


3 Analysis at 7 = —0.5 

3.1 Statistics 


We have used a vectorized three-hit Metropolis algorithm with an acceptance 
of about 50% at each hit. Due to the larger autocorrelation at 7 = —0.5 
compared with 7 = —0.2 the measurements have been done after 50 Metropolis 
sweeps each. 

Table 1 shows the lattice sizes L^ ■ L t with L t = 2 L s , the (3 values of the 
simulations, and the number of measurements in multiples of 1000. The phase 
transition turns out to be located at j3 c = 1.4223(2). 

The dagger f in table 1 denotes the simulations in which the two state signal 
is observed in the evolution of the plaquette energy. These runs have been 
excluded from the analysis. 
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Table 1 

Numbers of measurements in the confinement phase in multiples of 1000. All runs 
in which the two state signal was present are omitted from the anylsis. They are 
denoted by a f in this table. The phase transition is at (3 C = 1.4223(2). 


p 

16 3 32 

20 3 40 

24 3 48 

28 3 56 

32 3 64 

1.4 

4.5 





1.405 

5.0 





1.41 

3.6 





1.412 

3.0 





1.414 

3.6 





1.416 

3.0 





1.418 

4.1 

2.5 




1.419 

t 

5.5 

1.9 



1.42 

t 

5.7 

2.7 



1.4205 

t 

t 

4.1 

1.9 


1.4207 


t 

t 

t 

1.1 

1.4210 




t 

t 


The measurements of the GB masses use the techniques of smearing and diag- 
onalization of the correlation matrices described for pure QCD in [11,12], The 
measurement of the effective energies 6j(t) in the various GB channels and the 
way of estimating their statistical errors is explained in section 4 of [1], 


3.2 Gauge-ball channels 


The various gauge-ball channels result from the irreducible representations 
of the proper cubic group labeled R = A\, A 2 , E, T\, T 2 (with dimensions 
1,1, 2, 3, 3, respectively). Inversion then introduces a parity quantum num¬ 
ber P = ±1 and, since the representations are complex, one can also assign a 
charge parity C — ±1. The GB states are completely classified by R PC [13]. 

We measured the correlation functions of all the 20 channels R PC with mo¬ 
mentum p = 0 and some of them with the smallest lattice momentum p — 1 
(in units of 2tt/L s ). For the construction of operators with a certain symmetry 
see section 4 of [1] as well as [13,14], 

The effective energies €j (t) in each GB channel j = R P( have been determined 
for as many t < 6 as possible. The GB energies Ej can in principle be obtained 
from a fit to a plateau in the values of €j(t) which is expected for large t. 
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Fig. 1. Effective energies from A^ + (0) and T^^O) correlation functions (a) at 
/3 = 1.414 on the 16 3 32 lattice and (b) at j3 = 1.42 on the 24 3 48 lattice. The 
plateaus have been fitted omitting the points at t = 0,1 at /3 = 1.414, omitting 
t = 0,1 for Af + and t = 0,1,2 for T 2 ++ at j3 = 1.42. The resulting GB energies Ej 
and their errors are indicated by horizontal lines. 


Since the errors are sometimes very large or the effective energies could not 
be determined at all. the reliability of the determination of the energies Ej 
strongly depends on the GB channel. We classify this reliability in a similar 
way as done in [ 1 ]: 

+ + + states. At most /3 values a plateau in €j(t) is found and can be fitted 
at t > 2 or t > 3. 

++ states. At most (3 values a plateau in €j(t) is found only if t — 1 is 
included. 

+ states. Plateaus in 6j{t) can be found only at some f3 values and only if 
t = 1 is included. 

? states. No plateaus in 6j(t) are found: we cannot determine a GB energy 

at all. 

Table 2 shows the quality of the GB channels at 7 = —0.5. It is similar to 
the situation at 7 = — 0.2 [ 1 ] with some differences in the + and ? states. 
In contrast to [1] we have no (3 point close to the phase transition with a 
particularly high statistics. This prevents us from verifying small deviations 
from the one-particle dispersion relation for states denoted in [ 1 ] by an asterisk. 

The effective energies 6j(t) of Af + (+ + + quality) and T 2 ++ (++ quality) 
with p = 0 are shown in fig. 1 (a) at [3 = 1.414 on the 16 3 32 lattice and 
(b) at [3 = 1.42 (close to (3 C ) on the 24 3 48 lattice. The horizontal lines are the 
resulting GB energies Ej with errors. 

The GB masses rrij have been calculated from the GB energies Ej using the 
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Table 2 

GB states observed in various channels with symmetry R PC and momentum p = 
0,1 (in units of 2ir/L s ) and their scaling behaviour. If two momenta are listed in 
one line, the masses have been assumed to be the same. The ordering of the lines 
is the same as used in [1], The reliability of the energy determination is indicated 
as “quality”. In the columns “z/” and “mass” the scaling exponent of the GB mass 
and its approximate proportionality to one of the mass scales is indicated. The 
ratio of each mass to the corresponding mass scale is given in the column “ratio”. 
The ratios at 7 = —0.2 from [1] are shown for comparison. Since we used only GB 
masses mj(/3 ) < 1.0 in lattice units in the determination of the scaling behaviour, 
the amplitudes and therefore the ratios are not known for all of these channels. 


R PC {p ) 

Quality 

v Mass 

Ratio 

7 = -0.5 

Ratio 

7 = - 0.2 

^++( 0 , 1 ) 


i/ g m g 

1 

1 

Tf + (1) 



1.0(3) 

1 . 01 ( 6 ) 

^+-( 0 ) 



1 . 0 ( 1 ) 

1.00(4) 

E + ~( 0) 

+ + + 


1 . 0 ( 1 ) 

0.98(4) 

T+-(0,1) 


m ng 

1 

1 

T+-(0,1) 



1 . 0 ( 1 ) 

1.00(4) 

t 2 -( 1 ) 




1.04(5) 

E ++ ( 0,1) 

++ 

~ 2 m ng 

2 . 2 ( 2 ) 

2.06(9) 

T++(0,1) 



2 . 1 ( 2 ) 

2.07(9) 

^r + (o) 

? 



2.88(13) 

^r(o) 

? 

z/ ng 


2.73(14) 

^ + ( 0 ) 

+ 



3.22(16) 


+ 


3.0(3) 

2.73(15) 

E-( 0 ) 

+ 



2.95(14) 

T++(0,1) 

+ 


3.5(4) 

3.61(18) 

Tf (0,1) 

+ 

~ 3?n ng 

3.0(3) 

2.94(13) 

t 2 -+(o) 

+ 



3.10(14) 

t 2 - + ( 1 ) 

+ 



2.63(12) 

t 2 -(0) 

+ 



3.25(16) 

^+-( 0 ) 

+ 


3.0(3) 


^ 2 + ( 0 ) 

+ 




E~ + ( 0) 

+ 


2.8(3) 


Tf+(0) 

+ 


3.8(4) 
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lattice dispersion relation 


3 

2 (cosh. Ej — 1 ) = m* + 2 ^(1 — cos p^). ( 6 ) 

/•*= 1 


If the masses for a certain i? PC< are consistent for both p = 0 and p — 1, we 
assumed equality and used the same parameter in the fits. In this case both 
momenta are listed in the same line in table 2 and are denoted by the same 
symbol in figs. 2 and 4. 


3.3 Two mass scales 


Having determined the GB masses rrij at various (3, their scaling behaviour 
with f3 is analyzed. First, the power law 

m j = Cj(Pi j) ~ (3) Vl (7) 


is fitted to the individual GB masses measured in each GB channel j = R PC . 
The results for the critical values [3[p agree well within the error bars, so we 
assume equality [3p = /3 C and use the power law 

nrij = Cj((3 c — (3) Uj ( 8 ) 


with one /3 C for all GB channels. 

The resulting Uj from the fits ( 8 ) are shown in fig. 2. In these fits we restricted 
/3 to the values having rrij((3 ) <1.0 in lattice units. So the scaling behaviour 
could be determined for all GB channels with + + + and ++, but only for 
some of the + channels. 

Fig. 2 shows different critical exponents Vj for the scaling of GB masses in two 
groups of GB channels separated by the double vertical line in the figure. These 
two groups consist of the same GB channels as found at 7 = —0.2 [1], The 
two critical exponents are very similar to the results at 7 = —0.2, see (3),(4), 
and are therefore called u g and u ng respectively. The third column of table 2 
shows the association of GB channels to one of the two groups with u g or u ng . 

To demonstrate the difference between the exponents v g and z/ ng in a way 
independent of the numerical value /3 C , fig. 3 shows the masses of Af + and 
Ti~ in a doubly logarithmic plot. These two channels are the highest quality 
representatives of the u g and u Qg group, respectively. 
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Fig. 2. Values of Vj obtained separately in each channel. The double vertical line 
separates two groups with distinctly different vj . The dashed vertical lines separate 
channels with different quality within the group with the non-Gaussian exponent 
with decreasing quality from left to right. 
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7 

6 

5 


4 

2 3 4 5 6 7 

10 m(Aj + + ) 

Fig. 3. A log-log plot of the masses of the and gauge-balls in the interval 
1.412 < (3 < 1.4205. The results from the channels with p = 0,1 are denoted by 
circles and squares, respectively. The solid line is the result of a straight line fit 
considering errors in both directions. The dashed line with a gradient of one, i.e. 
assuming equal exponents, is shown for comparison. 
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Fig. 4. The amplitudes Cj of the fits (8) under the assumption that there are just 
two critical exponents. The vertical lines have the same meaning as in fig. 2. 

The slope of the solid line in fig. 3 is the ratio 

^ = 0.663(24). (9) 


The dashed line with slope one (if the exponents were equal) is shown for 
comparison. The result (9) is compatible with the value (5) obtained at 7 = 
- 0 . 2 . 

Assuming there are just two critical exponents u g and u ng we use the power 
law fit ( 8 ) again with only two different exponents Vj. Fig. 4 shows the results 
obtained for the amplitudes c r We find the same amplitude for all v g GB 
masses. The amplitudes of the v ng GB masses seem to be integer multiples of 
another amplitude value. This situation is again very similar to the results at 
7 = -0.2 [1], 

The ratios of the amplitudes Cj in each of the two groups of GB states have 
been determined using the amplitudes of the highest quality channels Af + 
and Tj +_ as a reference. These amplitude ratios determine the mass ratios in 
a possible continuum limit and are shown in the fifth column of table 2. The 
results at 7 = —0.2 from [1] are given for comparison. The amplitude ratios 
turn out to be compatible within error bars for all GB channels for which the 
amplitudes have been determined at both 7 values. 
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Fig. 5. Masses of the A^ + (lower points) and gauge-balls versus (3 for both 
momenta (p = 0,1 being denoted by circles and squares, respectively). The dashed 
curves are fits to the data shown. The full curves represent the mass scales obtained 
in a fit with all + + + GB states. 


To present results for the values of /3 C , u g , and u ng with reliable error estimates, 
we repeat the power law fit (8) with two exponents and fit the + + + gauge- 
balls only. The results are 


0.51(3), 

(10) 

0.35(2), 

(11) 


and 


I3 C = 1.4223(2). 


( 12 ) 


The ratio u ng / ,y g is consistent with (9). 

Fig. 5 shows the scaling of the masses of the GB states A± + and T 1 H with 
(3. The dashed curves are power law fits to the data shown while the solid 
curves represent the results from the fit to all + + + GB states which lead 
to (10)...(12). 
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.4 

.3 

.2 

1.400 1.405 1.410 1.415 1.420 

min 

Fig. 6 . Values of v g (circles) and v ng (squares) at 7 = —0.5 obtained in a fit with 
/dmin < (3 < (3 C (a) with (3 C = 1.4223 fixed and (b) with free (3 C . The horizontal lines 
indicate the values ( 10 ) and ( 11 ). 

4 Stability of v and comparison with essential singularity 


□ : u„ 




(a) 

y = -0.5 


.4 

.3 

.2 



1 400 1.405 1.410 1.415 1.420 

/^min 


4.1 Stability of v 


In order to test the reliability of the numerical determination of the exponents 
z/ g and z/ ng we repeat the power law fits ( 8 ) and restrict the (3 values to the 
range near the phase transition /3 min < (3 < (3 C . The resulting exponents are 
then plotted against /3 m i n to detect any possible non-asymptotic behaviour. 

Fig. 6 , part (a) shows the exponents u g and z/ ng against (3 m | n obtained from the 
power law fit ( 8 ) using the + + + GB masses with (3 C = 1.4223 fixed. In the 
fit with free (3 C in part (b) of fig. 6 we use only the masses in the two highest 
quality GB channels Af + and T+~ to avoid any influence caused by the fact 
that the masses of some GB states could be determined only at a few f3 values 
near the phase transition. 

The critical exponents v g) z/ ng , and even their statistical errors are very stable 
with respect to f3 m - m when (3 C is fixed (fig. 6 (a)). With free (3 C (fig. 6 (b)) the 
exponents are smaller if the fit is restricted to the [3 values near the phase 
transition, but their statistical errors grow by an amount comparable to this 
effect. 

For comparison we also present the same analysis at 7 = —0.2 using our data 
from [1], Fig. 7 (a) and (b) show these results with (3 C = 1.1609 fixed and with 
free (3 C using the same GB channels as in fig. 6 . 

The resulting critical exponents obtained with (3 C fixed turn out to be very 
stable. They are quite stable even with free f3 c . 

Currently we are unable to decide whether the decrease of the exponents at 
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Fig. 7. Values of n g (circles) and z/ ng (squares) at 7 = —0.2 obtained in a fit with 
/dmin < (3 < (3 C (a) with f3 c = 1.1609 fixed and (b) with free f3 c . The horizontal lines 
indicate the values (4) and (3). 



Fig. 8. Masses of the Af + (lower points) and gauge-balls (upper points) at 
7 = —0.2 (left plot) and 7 = —0.5 (right plot). The momenta p = 0,1 are denoted 
by circles and squares, respectively. The fits with the essential singularity law (13) 
with free /3 C using all data give the solid curves. The dashed curves are produced 
with fixed f3 c and restricting 1.152 < /3 < /3 C at 7 = — 0.2 resp. 1.416 < /3 < /3 C at 
7 = -0.5. 

7 = —0.5 with free (5 C near the transition is a statistically significant effect. 
The data at 7 = —0.2 have in general larger statistics than those at 7 = —0.5 
and do not show such a decrease (fig. 7 (b)). 

If the fits are made using g — g c instead of (3 — (3 C in ( 8 ), we find fully consistent 
values for the critical exponents z/ g and u ng . 


4-2 Comparison with essential singularity 

Motivated by [7] we have also explored the possible existence of an essential 
singularity and tried to describe the scaling of the masses of j = A± + and 
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j = T 1 + with the essential singularity law 


rrij(/3 ) = Cj 



Vj \ 


(13) 


with constants r]j. 

We studied several fits. Including all data and leaving all parameters of (13) 
free, the resulting /3 C is far off the anticipated value (solid curves in fig. 8 ). 
Fixing j3 c = 1.4223 (resp. f3 c = 1.1609 at 7 = —0.2) and restricting (13) to 
the data near the phase transition 1.416 < f3 < /3 C (resp. 1.152 < (3 < /3 C at 
7 = — 0 . 2 ) we obtain a satisfactory description of the masses (dashed curves). 
But in this case the results rjj are not very stable with respect to the choice of 
/3 C , e.g. varying /3 C = 1.4223 ± 0.0002 results in a variation of r/ g = 0.05 ± 0.02. 

It is thus possible that for data in the critical region the results are consistent 
with an essential singularity. However, establishing the existence of such a 
singularity would need much larger correlation lengths. 


5 Comparison of the results at 7 = —0.2 and 7 = —0.5 and conclu¬ 
sions 


We have compared the scaling of the gauge-ball masses when the phase tran¬ 
sition is approached at 7 = —0.5 from the confinement phase with the anal¬ 
ogous results obtained at 7 = —0.2 [1], Both the critical exponents and the 
amplitudes ratios of the GB masses which are expected to show a universal 
behaviour at a second order phase transition, have fully compatible values at 
both 7 . 

So it is most probable that the results obtained for the scaling of GB masses 
reflect properties of the line of deconfinement phase transitions in the param¬ 
eter space (/?, 7 ) of the extended Wilson action in a region of values 7 < —0.2. 
The most natural explanation of this universality would be the presence of 
a continuous phase transition in the coupling parameter space not far from, 
though not necessarily at the points we have investigated. Its localization may 
finally require the introduction of additional coupling parameters outside the 
(/3, 7 ) plane of the extended Wilson action. 

This is to be compared with other attempts to explain the occurrence of two 
distinct critical exponents. Certain scenarios (called T g and T ng in [1]) in which 
the cross-over near a tricritical point is responsible for the different exponents 
observed in non-vanishing distance from the critical line, are most probably 
ruled out as such an effect should be 7 -dependent. 
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The scenario suggested in [3], in which z/ ng is some effective value of v which 
decreases to 1/4 as the phase transition is approached, is not supported (but 
also not excluded) by our stability analysis. In any case z/ g ~ 0.5 is far away 
from the value 1/4 = 1/d expected for a first order transition. 

Of course, future work on larger lattices may disclose new facts. The question 
whether there is a certain range of parameters in which the phase transition is 
of second order is still open. However, our results obtained for the gauge-ball 
masses as small as 0.2 in lattice units are fully consistent with this possibility. 
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